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ABSTRACT 

We demonstrate that constraints on cosmological parameters from the distribution of clusters as a function of 
redshift (dN/dz) are complementary to accurate angular diameter distance (Da) measurements to clusters, and 
their combination significantly tightens constraints on the energy density content of the Universe. The number 
counts can be obtained from X-ray and/or SZ (Sunyaev-Zel'dovich effect) surveys, and the angular diameter 
distances can be determined from deep observations of the intra-cluster gas using their thermal bremsstrahlung 
X-ray emission and the SZ effect. We combine constraints from simulated cluster number counts expected from 
a 12 deg^ SZ cluster survey and constraints from simulated angular diameter distance measurements based on 
the X-ray /SZ method assuming a statistical accuracy of 10% in the angular diameter distance determination of 
100 clusters with redshifts less than 1.5. We find that $7,,, can be determined within about 25%, fl\ within 20%, 
and w within 16%. We show that combined dN/dz + D^ constraints can be used to constrain the different energy 
densities in the Universe even in the presence of a few percent redshift dependent systematic error in D^. We also 
address the question of how best to select clusters of galaxies for accurate diameter distance determinations. We 
show that the joint dN/dz+D^ constraints on cosmological parameters for a fixed target accuracy in the energy 
density parameters are optimized by selecting clusters with redshift upper cut-offs in the range 0.5 < z < 1. 

Subject headings: cosmological parameters - cosmology: theory - galaxiesxlusters: general 



1. INTRODUCTION 

A point often raised in the era of "precision cosmology" is 
that any particular method to determine fundamental cosmolog- 
ical parameters wiU suffer from degeneracies: only a combina- 
tion of parameters can be determined accurately (Bridle et al. 
2003). This issue has been highlighted recently by a number of 
papers focusing on the cosmological usefulness of galaxy clus- 
ters (Ettori, Tozzi & Rosati 2003; Vikhlinin et al. 2003; Mei & 
Bartlett 2003; Kujat et al. 2002; Newman et al. 2002; Rubino- 
Martin & Sunyaev 2002; Verde, Haiman & Spergel 2002; for 
a recent review based on X-ray observations see Rosati, Bor- 
gani & Norman 2002). Indeed, to date the tightest constraints 
have been derived using a combination of different methods, as 
demonstrated recently by the Wilkinson Microwave Anisotropy 
Probe (WMAP) team (Spergel et al. 2003). The importance of 
using independent methods to determine cosmological parame- 
ters is two-fold: (1) the methods may be combined to break de- 
generacies and better constrain individual parameters; and (2) 
consistency tests are possible and systematic errors can be stud- 
ied. 

Although we focus on constraints from the change in the 
number density of clusters as a function of redshift, dN/dz, and 
from measurements of the angular diameter distance of clusters 
based on the Sunyaev-Zel'dovich (SZ) effect. Da, further clus- 
ter observables could be used as additional indicators of cosmo- 
logical parameters. Levine et al. (2002) showed that using the 
cluster temperature function to constrain cosmological parame- 
ters can be made more efficient by adding the mass-temperature 
relation normalization from X-ray observations, and discussed 
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briefly how to optimize cluster surveys based on their method. 
Majumdar & Mohr (2003) emphasized that the mass-temperature 
relation with masses determined in a relatively modest follow- 
up program can be used in combination with the dN/dz test 
to greatly improve systematic limitations that arise from cluster 
evolution. Several other observables can be useful, such as scal- 
ing relations (Bialek et al. 2001 ; Verde et al. 2002), the shape of 
the cluster mass function (Hu 2003); and the three-dimensional 
cluster power spectrum (Refregier et al. 2002). 

Constraints on w as a function of redshift from planned SZ 
and supernovae surveys, and from their combination were stud- 
ied by Weller, Battye & Kneissl (2002) and WeUer & Albrecht 
(2002). Note, however, that cosmic variance can also be impor- 
tant in some cluster surveys (Hu & Kravtsov 2003). 

The redshift distribution of the number of clusters, dN/dz, 
is sensitive to both the change of the cosmological volume ele- 
ment and to the growth function of structure formation. Cluster 
number counts as a function of redshift can be determined from 
X-ray or SZ surveys (e.g. Carlstrom et al. 2002; Rosati et 
al. 2002). Haiman, Mohr & Holder (2001) discussed in detail 
constraints on dark energy (via the equation of state parameter, 
w) from the redshift distribution of clusters from future SZ and 
X-ray surveys and suggested that combining those constraints 
with constraints using Type la SNe or CMB fluctuations can hft 
degeneracies in the il,„-w plane (see also Wang & Steinhardt 
1998). Holder, Haiman & Mohr (2001) used the redshift dis- 
tribution of clusters from future SZ and X-ray surveys to con- 
strain il„j and r^A pointing out that, similarly to the parameter 
space (J7m,w), constraints from cluster redshift distribution are 
complementary to constraints from Type la SNe and CMB fluc- 
tuations. The X-ray thermal bremsstrahlung emission and SZ 
effect (Sunyaev & Zel'dovich 1980) depend on different combi- 
nations of the physical parameters of the cluster and cosmology 
and provide us with a way to determine the angular diameter 
distance. Da, to the cluster (the so-called X-ray/SZ method, 
for recent, detailed references, see Carlstrom, Holder & Reese 
2002 and Reese et al. 2003). The angular diameter distance 
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probes directly the curvature of the Universe. 

Molnar, Birkinshaw & Mushotzky (2002) discussed constraints 
on models with parameters {^lm,w,h), and {Q„,ilA,h) using sim- 
ulations of angular diameter distance measurements to clusters 
of galaxies and showed that the degeneracies in cosmological 
parameters from this technique are similar to those from SNe 
la, and therefore they are complementary to constraints from 
redshift distribution of clusters. This is not surprising, since 
the luminosity distance, which is utilized in the SN studies, is 
closely related to the angular diameter distance. 

In this paper, we show that it will be possible using clusters 
of galaxies alone to constrain accurately the matter density pa- 
rameter, rim, and the cosmological constant density parameter, 
J7a, or the equation-of-state parameter, w = p/p,hy combining 
dN/dz and Da measurements of clusters using the X-ray/SZ 
method, and estimate the accuracy achievable in these parame- 
ters. 

Since long observations are needed for high accuracy dis- 
tance measurements with present day observatories, a relatively 
small sample of clusters has to be selected. This cluster sam- 
ple does not have to be assembled from the same survey that is 
used for the cluster redshift distribution test. Any sufficiently 
accurate survey would suffice as long as it has the necessary 
redshift coverage. We address the natural question of how best 
to select samples of clusters of galaxies for efficient angular di- 
ameter distance determination. As an example, we discuss how 
we can optimize the selection of clusters for detailed distance 
measurements by allowing both the number of clusters in the 
sample and the upper redshift cut-off of the sample to be free 
parameters. We also discuss how to select a cluster sample to 
minimize the random errors. 

The rest of this paper is organized as follows. In § 2, we 
briefly summarize our methodology. In § 3, we apply our tech- 
nique to future cluster samples, and derive constraints on cos- 
mological parameters. In § 4, we discuss the issue of the op- 
timal selection of clusters for angular diameter distance mea- 
sures. Finally, in § 5, we offer our conclusions and surmnarize 
the implications of this work. 

2. METHODOLOGY 

We have studied combined constraints from the redshift dis- 
tribution of clusters and angular diameter distance measure- 
ments in two different sets of models. Both sets are described 
by the usual parameters (V,m,i^A,w,h,as,n), and we choose as 
our fiducial model (fim = 0.3,Qa = OJ,w = -l,h = 0.65, erg = 
0.9,n = 1). In Model A, we assume a cosmological constant 
(w = -1), but allow (il„,rt\) to vary independently. In Model B, 
we assume the universe is spatially flat (ilm + ^A = 1), but allow 
w to vary. We used constraints from previous simulations based 
on the redshift distribution of clusters derived from a deep 12 
deg^ SZ survey (see Holder et al. 2001, Model A; and Haiman 
et al. 2001, Model B). We then determined the expected con- 
straints on the same parameters from Da measurements with 
fixed accuracy, changing the total number of clusters and the 
upper redshift cutoff of the sample (for a detailed description of 
the method see Molnar et al. 2002). 

Most models of the X-ray emission from clusters of galaxies 
are based on the assumption of hydrostatic equilibrium. Note, 
however, that this is not a necessary assumption for the X- 
ray/SZ method to work. Detailed, spatially resolved X-ray 
spectroscopy can be used to observe the projected temperature 
and X-ray surface brightness distribution of the cluster and de- 



rive the de-projected temperature and gas density distribution, 
which, combined with the assumption of hydrostatic equilib- 
rium, can be used to determine the total mass distribution of 
the cluster (see for example Evrard et al. 2002). Recent XMM 
imaging spectra (Pratt & Arnaud 2003) support previous obser- 
vations of clusters which show that outside the central regions 
a large fraction of clusters are roughly isothermal and that the 
surface density can be described by a /3 model (Jones & For- 
man 1984; Mohr, Mathiesen & Evrard 1999), where the gas 
density can be written as p(r) = (1 + (r/rc)^)~^'^/^. The core ra- 
dius, rc, sets the scale of the cluster atmosphere and /3 speci- 
fies the slope of the density profile at large radius (Cavaliere & 
Fusco-Femiano 1976). For the sake of concreteness, we here 
assume an isothermal (3 model for the intra-cluster gas, how- 
ever this assumption is not essential in the determination of Da. 
The parameters that we need are: (1) and /3, which define 
the spatial structure of the intra-cluster gas and can be deduced 
from X-ray or SZ images; (2) the amplitude of the X-ray bright- 
ness and SZ decrement/increment (depending on the observing 
wavelength); (3) the temperature of the intra-cluster gas and 
the X-ray absorbing column density; and (4) the redshift of the 
cluster (derivable from x-ray or optical observations). 

The most important known redshift independent systematic 
errors in determining angular diameter distances arise from cal- 
ibration errors for the SZ and X-ray instruments, and finite 
cluster sizes from modeling the intra-cluster gas. Calibration 
errors cause systematic scale errors in the peak intensities of 
the SZ effect and X-ray surface brightness. The calibration er- 
rors of about 2.5% for interferometric observations (Carlstrom 
et al. 2002), and about 10% for XMM-Newton and Chandra, 
induce 5% and 10% errors in Da- The (3 model may give a di- 
vergent cluster mass if not truncated at some finite upper cut 
off radius. Assuming a j3 model with infinite extent introduces 
a systematic underestimate of Da- Based on numerical sim- 
ulations, Inagaki et al. (1995) estimated this systematic error 
to be as large as 10-20%. Systematic error estimates based 
on observations are much lower (up to 6%, Holzapfel et al. 
1997; Birkinshaw & Hughes 1994). However, in this paper 
we focus on the energy density and equation of state parame- 
ters which are not sensitive to redshift independent systematic 
errors since they depend only on the shape of the angular diam- 
eter distance-redshift function (for a detailed analysis see Mol- 
nar et al. 2002). 

Systematic errors depending on redshift do change the shape 
of the Da-z relation, and therefore could bias the inferred en- 
ergy density parameters. The most important such candidate is 
small scale gas clumping that varies with redshift. A general 
discussion of the effect of clumping on the determination of an- 
gular diameter distance can be found in Reese et al. (2003) and 
Molnar et al. (2002), and references therein. Numerical sim- 
ulations show that clumping caused by accretion shocks and 
mergers results in an overestimate in the angular diameter dis- 
tance to clusters by about 30% (LCDM, Mathiesen, Evrard & 
Mohr 1999). The error from clumping based on observations 
is estimated to be about 20% (Reese 2003). We will ignore 
this possibility here and simply note that (1) at present, there is 
no observational evidence for redshift dependent clumping, and 
(2) it is likely that a redshift dependent systematic effect would 
cause different systematic errors in dN/dz and Da measure- 
ments, and thus these kind of errors could be identified based 
on consistency. Overall a 30% error in the angular diameter 
distance is well within reach for a large sample of future sur- 
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Fig. 1 . — Simulated 97% confidence limits in the Q„, — Qf^ plane based on 
dN/dz and measurements. The outer set of solid, dashed, and dotted lines 
are constraints from the Da values alone for samples of clusters with upper 
redshift cutoffs of z = 1.5, 1.0, and 0.5. The dot-dashed line is the constraint 
from the dN/dz measured in a 12 deg^ SZ survey alone. The corresponding 
combined confidence intervals are shown as the irmer soUd, dashed, and dotted 
curves. 



veys (> 1000) using the spherical isothermal (3 model (Reese 
2003). We expect a much better accuracy for a smaller sample 
of clusters using spatially resolved spectroscopy. 

We assume further that the main known random error in de- 
termining Da, the orientation bias for non-spherical clusters, 
has been solved for our sample using galaxy velocity dispersion 
measurements, or a combination of X-ray, SZ, and weak lens- 
ing data with the assumption of hydrostatic equilibrium (Fox & 
Pen 2002). While this is an optimistic assumption, we find that 
in a sample of 100 clusters this effect would inflate our error 
estimates for and Oa by only an additional a = 5%. We 
expect about an 5% error in the cluster parameters from uncer- 
tainty in the spatial fitting, and a 5% error in the X-ray tem- 
perature is achievable with XMM-Newton and Chandra (Pratt 
& Arnaud 2003). For a detailed analysis of the error budget of 
the X-ray/SZ method, see Molnar et al. (2002), a sunnmary of 
errors based on observations can be found in Reese (2003). We 
assume an additional 5% statistical error from other sources. 
Finally we assume that the Hubble constant is determined from 
other measurements. 

Fox & Pen (2002), using numerical simulations, demonstrated 
that distances to clusters can be determined with a random error 
of about 5% combining X-ray, SZ, and lensing measurements 
and the assumption of hydrostatic equiUbrium. Although they 
ignored errors due to lensing, their result suggest that a 10% 
error in the determination of angular diameter distances is re- 
alistic in the near future. Therefore we performed simulations 
using a 10% error in the Da values. Note, however, that larger 
errors, as long as they are truly random, can be compensated 
with using more clusters, for example, assuming a 20% ran- 
dom error in Da, we would need a sample of about 400 clus- 
ters, and our final conclusion of this paper would be unaltered. 
In all of our simulations the redshift distribution was derived 
by randomly sampUng a uniform distribution in redshift space 
with an upper cut-off at fixed redshift. The number of sampled 
clusters in the redshift range 0.1-1.5 changes by less than a fac- 
tor of 3 as a function of cosmology, which, in practice, makes 
it easy to choose a quasi-uniform sample. Note, however, that 



Fig. 2. — Simulated 97% confidence limits in the Qm — w plane based on 
dN/dz and measurements. The line codes are the same as in Figure 1. 



our conclusions are not sensitive to the exact distribution of the 
number of clusters in redshift as long as they smoothly cover 
the targeted redshift range. 

Systematic effects in the cluster redshift distribution from the 
mass-X-ray temperature, the power spectrum normahzation, 
cTg, and the mass function were studied by Battye & Weller 
(2003). They concluded that these systematic effects are im- 
portant and more studies are necessary to reduce them. Sys- 
tematic effects in the cluster redshift distribution in SZ surveys 
from the surface brightness bias and cluster spatial orientation 
have not been discussed in detail. The detectibility is a function 
of signal/noise and, for extended sources (like clusters), this is 
determined by surface brightness and not by flux. X-ray and 
SZ surveys will not detect very extended clusters which have 
surface brightness too low to yield the necessary signal/noise. 
Even if we use a mass-size relation and a constant mass cut off, 
there is a dispersion around any observed or theoretical mass- 
size relation. 

As a result of the orientation bias, prolate clusters will be 
over-represented since they have larger central SZ and X-ray 
surface brightnesses than oblate clusters. Faint ellipsoidal clus- 
ters aligned with their long axis in the plane of the sky might 
not be detected, lower mass clusters with their long axis in the 
line of sight will be falsely detected in a given mass bin. How- 
ever, there are more low mass than high mass clusters, so we 
underestimate the number of clusters at the given mass thresh- 
old. Using a mass cut-off which is higher than the nominal 
detection-Umit-implied mass cut-off would give us a chance to 
estimate this bias. 

We used the likelihood function obtained by Holder et al. 
(2001) for model parameters ri„, and I^a (our Model A). Holder 
et al. defined their likelihood function based on the Cash C 
statistic (Cash 1979), kept h constant, and marginalized over 
the power spectrum normalization, o-g. We used the likelihood 
function of Haiman et al. (2001) for the cluster redshift distribu- 
tion test in parameter space and w (our Model B). Haiman 
et al. defined their likelihood function as the product of the 
Poisson probability of detecting a total number of clusters 
and a Kolmogorov-Smirnov probability for the unnormalized 
redshift distribution of these A'' clusters^. Haiman et al. also 

^ Note that using the likelihood function based on the Cash statistic would result 
in somewhat tighter constraints (as demonstrated expUcitly by a comparison of 
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Fig. 3. — Simulated 97% confidence limits in the f2„, — il^ plane based on 
dN jdz and Df^ measurements with redshift dependent systematic errors. The 
outer set of solid, dashed, and dotted lines are constraints from the Da values 
alone for a sample of 100 clusters with random errors of 10%, and redshift 
dependent systematic errors of +5% and -5% at z=l growing linearly from 
at z=0. The dot-dashed line is the constraint from the dN jdz measured in a 
12 deg^ SZ survey alone. The corresponding combined confidence intervals 
are shown as the iimer solid, dashed, and dotted curves. 



kept h constant and marginalized over the power spectrum nor- 
malization. We kept h constant, assumed Gaussian errors in the 
angular diameter distance measurements, and constructed the 
hkehhood function based on the usual Ax^ statistic. The an- 
gular diameter distance function does not depend on the power 
spectrum, so we can combine these two likelihood functions in 
a straightforward way by multiplying them. In more detailed 
work on specific observational strategies it would be necessary 
to multiply the hkehhood functions before marginalizing over 
unconstrained parameters. 

3. COMBINED CONSTRAINTS ON COSMOLOGICAL 
PARAMETERS 

In our first set of simulations, we assumed 100 clusters with 
an error of 10% in the Da determinations with different redshift 
cut-offs out to z= 1.5. We show the 97% confidence hmits 
(CLs) for the combined constraints from the redshift distribu- 
tion of clusters and angular diameter distance measurements 
based on simulations in Figures 1 and 2. In Figure 1 we show 
constraints on fi^ and f^A (Model A). This figure shows that, 
as we select clusters with higher and higher redshift cut-offs, 
the constraints from D,\ measurements alone get tighter and the 
confidence contours rotate counter-clockwise. The constraints 
get tighter because at higher redshifts, T>a(j) is more sensitive to 
these parameters, and rotate, since these parameters are sensi- 
tive to different combinations of the density parameters at dif- 
ferent redshifts. As a consequence of this rotation, although 
at higher upper redshift cut-offs the constraints from Da mea- 
surements are tighter, they are farther from orthogonal to con- 
straints from the redshift distribution of clusters. The combined 
constraints do not improve substantially from upper cut-offs of 
about z> I , although the likelihood function does get somewhat 
narrower. 

We conclude that, contrary to naive expectations, observa- 
tions of clusters at redshifts exceeding z » 1 wiU not improve 
the constraints on and using these methods. Rather, the 

these two different statistics in Verde et al. 2001). 
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Fig. 4. — Simulated 97% confidence limits in the Qm — w plane based on 
dN /dz and measurements with redshift dependent systematic errors. The 
line codes are the same as in Figure 3. 



most important parameter is simply the total number of clusters, 

as long as we observe clusters with redshifts up to z '--^ 1 . This is 
good news because, while it is difficult to determine distances 
at high redshifts (z^l) due to the low flux of these clusters in 
the X-ray band, recent XMM-Newton and Chandra observa- 
tions readily find clusters out to redshifts of about z = 0-9 (see 
for example Ettori et al. 2002). Based on our simulations we 
find that the combined constraints from 100 clusters uniformly 
distributed with an upper cut-off at z = 1 could be used to de- 
termine Qm to within about 25%, and Oa to within 20% (97% 
CL). 

We show our results on combined constraints on 0,„ and w 
(Model B) from simulations in Figure 2. The figure shows that 
the constraints from angular diameter distances alone are once 
again getting tighter as we choose higher redshift cut-offs, but 
the constraints also spread more in w (they cross the w axis at 
higher values of w). As a result, as it can be seen from this 
figure, the high-w hmit does not improve by choosing clusters 
with redshift cut-offs greater than about z ~ 1 . Using 100 clus- 
ters, the combined method could determine w to about 16% 
(97% CL). 

We choose to study constraints on the energy densities in the 
Universe because they are not sensitive to systematic effects 
which are redshift independent. However, the energy densi- 
ties are sensitive to redshift dependent systematics. Although 
such an effect has not yet been found, a few percent redshift 
dependent systematic error would be difficult to identify obser- 
vationaUy. We carried out simulations to study how constraints 
on Ct„ and (Model A) and on 0,„ and w (Model B) change 
assuming redshift dependent systematics at few percent level. 
We show our results on Figures 3 and 4 assuming a systematic 
bias in the inferred distance that grows linearly from no bias at 
z=0 to a +(-)5% over(under)estimate of all distances at z=l us- 
ing 100 clusters and 10% random errors, as before. Combining 
the Da and dN /dz results still permits the measurement of the 
densities of matter and dark energy, or the density of matter and 
the equation of state parameter. Although a systematic shift in 
the derived values of these parameters is caused, such 5% sys- 
tematic errors do not move the error ellipse away from the in- 
put model values significantly. While larger redshift-dependent 
systematic errors of this type would lead to significant system- 
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atic errors in the derived parameters, they would also cause in- 
consistencies in the Da and dN/dz tests through which they 
could be detected. 

4. OPTIMIZING THE SAMPLE SELECTION FOR DISTANCE 

MEASUREMENTS 

Although we expect hundreds of clusters to be discovered in 
future X-ray and SZ surveys (for example, the XMM-Newton 
64 deg^ Large Scale Structure Survey will provide hundreds 
of candidate clusters, Refregier et al. 2002), it is not feasible 
to determine accurate distances to all discovered clusters. The 
reason for this is two-fold: (1) the observational time necessary 
for accurate distance determination rapidly increases with dis- 
tance, and is on the order of 50-100 ks with XMM-Newton and 
Chandra at redshifts of about 0.7; and (2) not all clusters are 
relaxed. 

The question then naturally arises: How can we assemble our 
cluster sample optimally, yielding the most precise and robust 
constraints on cosmological parameters? A fundamental physi- 
cal selection criterion is that the clusters should be as relaxed as 
possible. Based on numerical simulations, Roettiger, Stone & 
Mushotzky (1997), in order to minimize the effects of dynam- 
ical activity, suggest to exclude clusters with non-cylindrical 
surface brightness and temperature distribution in projection, 
with twisted isophots, with anisotropic galaxy velocity distri- 
bution, and those clusters for which /3 from spatial fitting is sig- 
nificantly different from (5 determined from spectroscopy. If we 
wish to eliminate the orientation bias, we need to use clusters 
which are, in most part, in hydrostatic equilibrium. Although 
the X-ray/SZ method works for any cluster as long as we have 
an accurate physical model for the gas, clusters which are not in 
dynamical equilibrium are difficult to model, contribute a large 
systematic error, and so should be excluded. The signature of 
such objects is complicated morphology of the spatial structure 
in the X-ray surface brightness and/or temperature image. Line 
of sight mergers would be difficult to identify based on mor- 
phology. Mergers might be identified by comparing distances 
determined using the X-ray/SZ and other methods. The veloc- 
ity distribution of member galaxies can also be used to identify 
fine of sight merging and numerical simulations can be used to 
verify merging (Gomez, Hughes & Birkinshaw 2000). A clus- 
ter which is an outlier in scaling relations (X-ray temperature - 
luminosity, mass - X-ray temperature, mass - X-ray luminos- 
ity) would be a likely case of dynamical activity, since on-going 
merging leads to temporary enhancements of the X-ray flux and 
raised average intra-cluster gas temperatures because of shock 
heating (Randall, Sarazin & Ricker 2002). 

In our previous simulations, we kept the number of clusters 
fixed at N = 100, in order to see whether constraints improve 
by using clusters at z > 1. Since it takes longer to achieve the 
target accuracy at higher redshifts, one has to consider the ob- 
servational time as well. A detailed analysis of how the required 
observational time for a given Da accuracy depends on the red- 
shift of the cluster is beyond the scope of this paper since the 
necessary parameters are instrument specific. Instead we in- 
vert the question, and use simulations to determine the number 
of clusters necessary to observe to achieve a fixed accuracy on 
the cosmological parameters, as a function of the upper redshift 
cut-off for the Da sample. For concreteness, we fix the errors at 
the level achievable with observations of 100 clusters with 10% 
accuracy in the angular diameter distance determination with 
an upper cut-off of z = 1, i.e. we determine the number of clus- 
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Fig. 5. — The number of clusters, N, in the sample for angular diameter 
distance determination that is necessary to achieve a fixed statistical accuracy, 
25% in 20% in Ua (Model A) and 16% in w (Model B) from combining 
SZ survey and Z)^ measurements, as a function of the upper redshift cut-off 
of the Da sample. The squares and the solid curve represent the number of 
clusters necessary to determine Clm and CI a, the triangles and the dashed curve 
represent the number of clusters needed to determine w. 



ters necessary to achieve an accuracy of about 25% in $1^, and 
20% in r^A (for our Model A), and 16% in w (for our Model B) 
as a function of upper redshift cut-off in the Da cluster sample. 

In general, different numbers of clusters are necessary to 
achieve any specified accuracy for different parameters . In prac- 
tice, we find that the ratio of errors does not change significantly 
in the considered redshift range and therefore once the ratio be- 
tween errors in fl^ and flA is fixed, it is possible to determine 
the number of clusters in the sample to achieve the same ratio 
in accuracy in both parameters as a function of redshift cut- 
off. Since adding constraints from Da measurements to those 
from the redshift distribution of clusters in Model B does not 
improve constraints on the matter density, we determine only 
the number of clusters necessary to achieve a fixed accuracy in 
w as a function of redshift. We show our results in Figure 5. In- 
terestingly, the curves give rather similar requirements for the 
number of clusters. The figure also clearly shows that, as ex- 
pected, with low upper cut-offs on redshifts of the sample, the 
necessary number of clusters increases, since the angular diam- 
eter distance is insensitive to the energy content of the Universe 
at low redshift. In accordance with our conclusions from Fig- 
ures 1 and 2, Figure 5 shows that at redshift cutoffs z>l the 
number of clusters needed tends to a fixed value. However, the 
observing time will increase significantly if one wants to main- 
tain the assumed 10% accuracy in the angular diameter distance 
measurements. Overall, Figure 5 suggests that the optimal red- 
shift range of a cluster sample for useful Da measurements is 
0.5<z<l. 

5. CONCLUSIONS 

We have demonstrated that the energy content of the Uni- 
verse can be constrained to high statistical accuracy using the 
cluster redshift distribution and angular diameter distance mea- 
surements. Degeneracies in cosmological parameters that are 
constrained by either observable are substantially weakened when 
they are used in tandem. We quantified constraints from the 
simulated cluster redshift distribution expected from a 12 deg^ 
SZ cluster survey and constraints from simulated angular di- 
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ameter distance measurements based on using the X-ray/SZ 
method, assuming an expected accuracy of 10% in the angu- 
lar diameter distance determination of 100 clusters with red- 
shifts z<l .5. We find that flm can be determined to a statistical 
accuracy (97% CL) of about 25%, within 20%, and w to 
an accuracy of about 16%. We also addressed the question of 
how to select clusters of galaxies optimally for accurate diam- 
eter distance determinations. Our results indicate that the joint 
dN/dz+DA constraints on cosmological parameters for a given 
observation time are optimized by using cluster Da measure- 
ments in the redshift range 0.5<z<l. We carried out simula- 
tions to study how combined constraints on the energy densities 
depend on redshift dependent systematic errors. We found that 
a combination of cluster redshift distribution and angular di- 
ameter distance determination measurements with an assumed 
redshift dependent systematic error which grows linearly to 5% 
at z = 1 still leads to significantly improved constraints on cos- 
mological parameters, without unacceptable systematic errors 
in those parameters. 

Comparison of the errors derived from the dN/dz+DA tech- 
nique with 100 clusters compares favorably with those from 
simulations based on type la supemovae (e.g. Gerke & Efs- 
tathiou 2002; see also Hannestad & Mortsell 2002; Perlmut- 
ter et al. 1999), with similar errors on w, but significantly 
tighter errors on Inclusion of priors from the WMAP ex- 
periment will tighten confidence limits further (complementar- 
ity between cosmic microwave background and cluster dN/dz 
constraints have been explicitly demonstrated in Haiman et al. 
2001). We beheve that the results in this paper are important, 
because they indicate that independent, tight statistical con- 
straints on cosmological parameters will be available internally 
from galaxy cluster surveys. Having independent constraints 
from different methods will be invaluable in understanding the 
systematic errors that will Ukely hmit our determination of cos- 
mological parameters in the next several years. 
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NASA LTSA Grant NAG5-3432 (PI: J. R Hughes). 
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